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ABSTRACT 

Severe acute respiratory syndrome coronavirus (SARS-CoV) causes a respiratory disease with a mortality rate of 10%. A mouse- 
adapted SARS-CoV (SARS-CoV-MA15) lacking the envelope (E) protein (rSARS-CoV-MA15-AE) is attenuated in vivo. To iden- 
tify E protein regions and host responses that contribute to rsARS-CoV-MA15-AE attenuation, several mutants (rSARS-CoV- 
MA15-E*) containing point mutations or deletions in the amino-terminal or the carboxy-terminal regions of the E protein were 
generated. Amino acid substitutions in the amino terminus, or deletion of regions in the internal carboxy-terminal region of E 
protein, led to virus attenuation. Attenuated viruses induced minimal lung injury, diminished limited neutrophil influx, and 
increased CD4* and CD8* T cell counts in the lungs of BALB/c mice, compared to mice infected with the wild-type virus. To 
analyze the host responses leading to rSARS-CoV-MA15-E* attenuation, differences in gene expression elicited by the native and 
mutant viruses in the lungs of infected mice were determined. Expression levels of a large number of proinflammatory cytokines 
associated with lung injury were reduced in the lungs of rsARS-CoV-MA15-E*-infected mice, whereas the levels of anti-inflam- 
matory cytokines were increased, both at the mRNA and protein levels. These results suggested that the reduction in lung in- 
flammation together with a more robust antiviral T cell response contributed to rsARS-CoV-MA15-E* attenuation. The attenu- 


ated viruses completely protected mice against challenge with the lethal parental virus, indicating that these viruses are 


promising vaccine candidates. 


IMPORTANCE 


Human coronaviruses are important zoonotic pathogens. SARS-CoV caused a worldwide epidemic infecting more than 8,000 
people with a mortality of around 10%. Therefore, understanding the virulence mechanisms of this pathogen and developing 
efficacious vaccines are of high importance to prevent epidemics from this and other human coronaviruses. Previously, we dem- 
onstrated that a SARS-CoV lacking the E protein was attenuated in vivo. Here, we show that small deletions and modifications 
within the E protein led to virus attenuation, manifested by minimal lung injury, limited neutrophil influx to the lungs, reduced 
expression of proinflammatory cytokines, increased anti-inflammatory cytokine levels, and enhanced CD4* and CD8* T cell 
counts in vivo, suggesting that these phenomena contribute to virus attenuation. The attenuated mutants fully protected mice 
from challenge with virulent virus. These studies show that mutations in the E protein are not well tolerated and indicate that 


this protein is an excellent target for vaccine development. 


bubs (CoVs) are responsible for a wide range of im- 
portant veterinary and human diseases (1). Severe acute res- 
piratory syndrome (SARS) is caused by a coronavirus (SARS- 
CoV) that emerged in Guangdong Province, China, causing the 
2002-2003 epidemic, and infecting more than 8,000 individuals 
with a 10% mortality rate (1-5). A novel human coronavirus, 
Middle East respiratory syndrome coronavirus (MERS-CoV), was 
isolated from the sputum of a man with pneumonia and renal 
failure in Jeddah, Saudi Arabia, in 2012 (6, 7). As of 21 November 
2014, the WHO confirmed 909 MERS cases, in which 36% of 
patients died (8). CoVs similar to SARS-CoV and MERS-CoV 
have been found in bats distributed across the world from which 
new zoonotic transmissions into human population could occur 
(9-15). Therefore, understanding the virulence mechanisms of 
these pathogens and developing efficacious vaccines and therapies 
are of high importance. 

SARS-CoV is an enveloped virus with a single-stranded, posi- 
tive-sense, 29.7-kb RNA genome (16, 17). The 5’ two-thirds of the 
genome comprise two overlapping open reading frames, ORF la 
and ORF 1b, encoding two polyproteins, ppla and pp1lab. In ad- 
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dition, SARS-CoV encodes a set of structural proteins present in 

all CoVs: spike (S), membrane (M), envelope (E), nucleocapsid 

(N), and accessory proteins 3a, 6, 7a, 7b, 8a, 8b, and 9b (18). 
SARS-CoV E protein is a small integral membrane protein of 
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TABLE 1 Primers used for the generation of rsARS-CoV-MA15-E* protein deletion mutants 


E* mutant PCRno. Primer Sequence (5’—3') 
A2 mutant, 24nt* 1 E-SARS-25871-VS = CGTTGTACATGGCTATTTCACCG 
1 E-SARS-26278-RS —§ GGTTTTACTAAACTCACGTTAACAATAAGCGCAGTAAGGATGGCTAGTGTG 
2 E-SARS-26223-VS § GCGCTTATTGTTAACGTGAGTTTAGTAAAACC 
2 E-SARS-28160-RS § CTGAGTGAGCTGTGAACC 
A3 mutant, 21 nt 1 E-SARS-25871-VS = CGTTGTACATGGCTATTTCACCG 
1 E-SARS-26299-RS —§ CGCGAGTAGACGTAAACCGTTGGTTTATTGCAGCAGTACGCACACAATCG 
2 E-SARS-26274-VS § AAACCAACGGTTTACGTCTACTCGCG 
2 E-SARS-28160-RS |= CTGAGTGAGCTGTGAACC 
A4 mutant, 18 nt 1 E-SARS-25871-VS § CGTTGTACATGGCTATTTCACCG 
1 E-SARS-26298-RS —§ CGCGAGTAGACGTAAACCACGTTAACAATATTGCAGCAGTACGC 
2 E-SARS-26248-VS § GCGTACTGCTGCAATATTGTTAACGTGGTTTACGTCTACTCGCG 
2 E-SARS-28160-RS © CTGAGTGAGCTGTGAACC 
A5 mutant, 36 nt 1 E-SARS-25871-VS § CGTTGTACATGGCTATTTCACCG 
1 E-SARS-26332-RS —§ GGAACTCCTTCAGAAGAGTTCAGTACTAAACTCACGTTAACAATATTGC 
2 E-SARS-26266-VS =GTTTAGTACTGAACTCTTCTGAAGGAGTTCC 
2 E-SARS-28160-RS © CTGAGTGAGCTGTGAACC 
A6 mutant, 36 nt 1 E-SARS-25871-VS = CGTTGTACATGGCTATTTCACCG 
1 E-SARS-26379-RS © CCAAACAGAATAATAATAATAGTTAGTTCGTTTAATTTTTAACACGCGAGTAGACGTAAACCG 
2 E-SARS-26296-VS © CGCGTGTTAAAAATTAAACGAACTAACTATTATTATTATTCTGTTTGG 
2 E-SARS-28160-RS |» CTGAGTGAGCTGTGAACC 


“ nt, nucleotides. 


76 amino acids. This protein contains a short hydrophilic amino 
terminus, a hydrophobic region, and a hydrophilic carboxy ter- 
minus (19). The hydrophobic region forms one amphipathic 
a-helix that oligomerizes and displays ion channel activity (20- 
22). E protein is present in small amounts within viral particles but 
is abundantly synthesized in infected cells (23-25), where it local- 
izes mainly in the endoplasmic reticulum Golgi intermediate 
compartment (ERGIC), participating in virus morphogenesis and 
budding (25-27). SARS-CoV lacking the E gene (rSARS-CoV- 
AE) was attenuated in vivo (28, 29) and protected against chal- 
lenge with virulent SARS-CoV (30-32). The E protein is a viru- 
lence factor, regulating cell stress response and apoptosis and 
promoting inflammation (33). 

SARS-CoV causes a respiratory illness characterized by acute 
lung injury (ALI), and its most severe pathological form is acute 
respiratory distress syndrome (ARDS). These pathological condi- 
tions are characterized by diffuse alveolar damage (DAD), pulmo- 
nary cellular infiltration, hyaline membrane formation, and 
edema accumulation, leading to hypoxemia and eventually to 
death (34-37). SARS-CoV infection induces migration to the 
lungs of immune cells, such as neutrophils and macrophages, 
which initiate and amplify the inflammatory response (38, 39). 
Accordingly, the levels of several proinflammatory cytokines, in- 
cluding gamma interferon (IFN-y), IFN-inducible protein 10 (IP- 
10)/CXCL10, monocyte chemoattractant protein 1 (MCP-1)/ 
CCL2, and several interleukins (IL-1, IL-6, IL-8, and IL-12), are 
elevated in the lungs and peripheral blood of SARS patients (40— 
49). This exacerbated inflammatory process correlates with lung 
injury and a poor outcome. T cells are essential to resolve SARS- 
CoV infections. T cell responses play a crucial role in SARS-CoV 
clearance and in protection from clinical disease (50-52). Accord- 
ingly, one notable finding in human SARS, associated with an 
adverse outcome, was the rapid development of lymphopenia, 
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with numbers of CD4* T cells more severely reduced than those of 
CD8* T cells during acute disease (38, 53-56). 

The SARS-CoV Urbani strain virus causes no significant dis- 
ease in wild-type (wt) mice (57). Passage through BALB/c mouse 
lungs resulted in a mouse-adapted virus (MA15 strain) (57), 
which upon infection reproduced many aspects of the human 
disease, such as high virus titers, pathological changes in lungs, 
viremia, neutrophilia, and lethality (57). We previously showed 
that rsARS-CoV-MA15-AE was attenuated in mice and that im- 
munization with rsARS-CoV-MA15-AE completely protected 
young and aged BALB/c mice against challenge with a lethal dose 
of MA15 (32, 58). Furthermore, we showed that SARS-CoV E 


TABLE 2 TaqMan assays used to analyze the expression of cellular 
genes by quantitative RT-PCR 


Proteinname TaqMan assay* Description 

TNE Mm00443258-ml Tumor necrosis factor 
IL-4 Mm00445259 ml Interleukin 4 

IL-5 Mm00439646_ml Interleukin 5 

IL-6 Mm00446190-ml Interleukin 6 

IL-10 Mm00439614 ml Interleukin 10 

IL-12B Mm00434174_ml Interleukin 12B 

IL-13 Mm00434204 ml Interleukin 13 


CXCL1/NAP-3 
CXCL2/MIP-2 
CXCL10/IP-10 


Mm04207460-ml 
Mm00436450-ml 
Mm00445235-ml 


Neutrophil-activating protein 3 
Macrophage inflammatory protein 2 
Interferon-inducible protein 10 


CCL2/MCP-1 Mm00441242-ml Monocyte chemotactic protein 1 
CCL3/MIPIA Mm00441259_gl Macrophage inflammatory protein la 
CCL4/MIP1B =Mm00443111_ml Macrophage inflammatory protein 18 
IEN-y Mm01168134-ml Gamma interferon 

TGF-8 Mm01178820_ml Transforming growth factor beta 1 
18S Mm03928990-gl_ — Ribosomic RNA 18S 


@Mm, Mus musculus. 
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FIG1 Schematic of mutations and deletions introduced within SARS-CoV E protein and growth kinetics of the mutant viruses (rSARS-CoV-MA15-E*). (A) The 
SARS-CoV genome is shown in the top, and the expanded region shows the E protein sequence and its different regions. White boxes represent the amino acids 
deleted within the E protein in each virus. Gray letters indicate the amino acids mutated to change the amino-terminal region of the protein. (B) Mutant virus 
growth kinetics. Subconfluent monolayers of Vero E6 and Huh7.5.1 cells were infected with wt, AE, and rsARS-CoV-MA15-E* viruses at an MOI of 0.001 on 
Vero E6 cells. Error bars represent standard deviations of the mean using data from three independent experiments. Hexagon, WT; diamond, AE; pentagon, Mut 


1; triangle, A2; star, A3; inverted triangle, A4; square, A5; and circle, A6. 


protein ion channel activity promoted virus virulence and fitness 
(37). To identify the E protein regions and the mechanisms lead- 
ing to SARS-CoV-AE attenuation, several mouse-adapted virus 
mutants encoding amino acid substitutions in the amino terminal 
or small deletions located in the carboxy-terminal region of the E 
protein (rSARS-CoV-MA15-E*) were constructed. Amino acid 
substitutions in the amino terminal, or deletion of regions in the 
central carboxy-terminal region of E protein, resulted in virus 
attenuation, accompanied by reductions in lung inflammation, 
neutrophil influx into the lungs, and proinflammatory cytokine 
expression. Remarkably, the number of T cells was increased in 
the lungs of mice infected with the less pathogenic viruses, most 
probably contributing to their more rapid clearance. Importantly, 
the attenuated mutants protected against the challenge with the 
virulent wt virus. 
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MATERIALS AND METHODS 


Cells. Vero E6, BHK, and Huh7.5.1 cells were kindly provided by E. Sni- 
jder (University of Leiden, The Netherlands), H. Laude (Unité de Virolo- 
gie et Immunologie Molecularies, INRA, France), and R. Bartenschlager 
(Department of Molecular Virology, University of Heidelberg, Germany), 
respectively, and were propagated as described previously (28). 

Viruses. Mouse-adapted SARS-CoV-MA15 (57) wasa gift from Kanta 
Subbarao (National Institutes of Health, Bethesda, MD). Recombinant 
viruses were rescued from infectious cDNA clones generated in our labo- 
ratory (32, 58, 59). 

Mice. Specific-pathogen-free BALB/c mice were purchased from the 
National Cancer Institute at the age of 6 or 16 weeks or from Harlan 
Laboratories (Holland) at the age of 8 weeks and maintained for 8 addi- 
tional weeks. All experiments involving SARS-CoV were conducted in 
biosafety level 3 laboratories in the animal care facility at the University of 
Iowa or at the Center for Animal Health Research (CISA-INIA, Spain), 
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FIG 2 Subcellular localization of mutant E proteins. (A) Vero E6 cells were infected with either rsARS-CoV-MA15-AE, -A3, or -A5 or wt recombinant viruses, 
at an MOI of 0.3, and fixed at 24 hpi. E protein (green) and ERGIC (red) were labeled with specific antibodies. Nuclei were stained with DAPI (blue). Merge 
indicates superposition of both labels. Original magnification was X 126. (B) The panel represents the percentage of overlap coefficient between E protein and 


ERGIC, calculated with Leica LAS AF v2.6.0 software. 


equipped with ventilated racks (animal transport unit and biocontain- 
ment unit; Allentown, Inc.) to store the animals during the experiment. 
All the protocols were approved by the EU, by the CISA-INIA Com- 
mittees of Animal Care Biosecurity and Bioethics, or by the University 
of Iowa Animal Use and Care Committee. All personnel were equipped 
with positive-pressure air-purifying respirators (3M HEPA AirMate, 
St. Paul, MN), 

Construction of pBAC-SARS-CoV-MA15-E* plasmids. Mutant vi- 
ruses (rSARS-CoV-MA15-E*) with amino acid substitutions in the amino- 
terminal region (rSARS-CoV-MA15-Mut 1) or with small deletions 
covering different regions of the carboxy terminus of the E protein 
(rSARS-CoV-MA15-A2, -A3, -A4, -A5, and -A6) were constructed using 
an infectious cDNA clone. cDNA encoding the genome of the SARS-CoV- 
MA\IS5 strain was assembled in a bacterial artificial chromosome (BAC) 
(pBAC-SARS-CoV-MAI15 plasmid) (32, 58). DNA fragments containing 


L REP 1a 


REP 1b 


nucleotides (nt) 26044 to 26779 of the SARS-CoV genome were generated 
by overlap extension PCR in the case of the carboxy-terminal mutants, 
using as the template the pBAC-SARS-CoV-MA15 plasmid and the prim- 
ers indicated in Table 1, and by gene synthesis (Bio Basic, Inc.) in the case 
of the amino-terminal mutant. This last fragment included four point 
mutations, generating four amino acid changes: S3A (TCA to GCA), V5L 
(GTT to CTT), T9A (ACA to GCA), and T11A (ACG to GCG). The final 
PCR products or the fragments generated by gene synthesis were digested 
with the enzymes BamHI and Mfel and cloned into the intermediate 
plasmid psl1190+BamHI/SacII SARS-CoV to generate plasmids 
psl1190+BamHI/SacII SARS-CoV-E* (psl1190-Mut 1, A2, A3, A4, A5, 
and A6). Plasmid psl1190+ BamHI/SaclII SARS-CoV contains a fragment 
corresponding to nucleotides 26045 to 30091 of the SARS-CoV infectious 
cDNA clone (59) engineered into the psl1190 plasmid (Pharmacia) using 
unique BamHI and SaclI restriction sites. Finally, these fragments were 
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FIG 3 Stability of SARS-CoV-MA15-E* after serial infections. The stability of rSARS-CoV-MA15-E* virus deletion mutants was examined after 8 passages in 
Vero E6 cells by sequence analysis. Asterisks denote the presence of nucleotide substitutions. 
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FIG 4 Weight loss and survival rate of mice inoculated with rsARS-CoV-MA15-E* mutants. BALB/c mice were intranasally infected with 1 X 10° PFU of each 
virus (n = 5 mice). Animals were monitored daily for weight loss (A) and survival (B). Animals that lost more than 30% of their initial body weight were 
euthanized. Differences in weight loss between attenuated and virulent viruses were statistically significant (*, P< 0.05). 


inserted into pBAC-SARS-CoV-MAI15 to generate pBAC-SARS-CoV- 
MA15-E* plasmids. The viruses were rescued in BHK and Vero E6 cells as 
previously described (28). Virus was cloned by three rounds of plaque 
purification. 

rSARS-CoV-MA15-E* growth kinetics. Subconfluent monolayers 
(90% confluence) of Vero E6 or Huh7.5.1 cells were infected at a multi- 
plicity of infection (MOI) of 0.001 with rsARS-CoV-MA15-E*, rSARS- 
CoV-MAI15-AE, or rsARS-CoV-MA15. Culture supernatants were col- 
lected at the indicated times postinfection, and titers were determined on 
Vero E6 cells as previously described (28). 

Indirect immunofluorescence microscopy. To detect viral E protein 
expression, Vero E6 cells were grown to 70% confluence on glass cover- 
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slips and infected with rsARS-CoV-MA15, rSARS-CoV-MA15-AE, and 
rSARS-CoV-MA15-E* at an MOI of 0.3. At 24 h postinfection (hpi), cells 
were fixed with 4% paraformaldehyde for 20 min at room temperature. 
After washing and permeabilization, primary antibody incubations were 
performed for 90 min at room temperature. Rabbit polyclonal antibody 
(Ab) specific for E protein (dilution of 1:500; kindly provided by Shen 
Shuo, Institute of Molecular and Cellular Biology, Singapore) and a 
mouse anti-ERGIC marker (dilution of 1:200; Alexis Biochemicals) were 
used. Secondary antibodies specific for rabbit or mice species conjugated 
to Alexa 488 (for detection of E protein) or Alexa 594 (for detection of 
ERGIC), respectively (dilution of 1:500; Invitrogen), were incubated for 
45 min at room temperature. Nuclei were stained using 4’ ,6-diamidino- 
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2-phenylindole (DAPI) (dilution of 1:200; Sigma), and slides were exam- 
ined on a Leica SP5 confocal microscope (Leica Microsystems). 

Specific infectivity of rsARS-CoV-MA15-E* mutants. Subconfluent 
monolayers (90% confluence) of Vero E6 cells were infected at an MOI of 
0.3 with rsARS-CoV-MA15-E*, rsARS-CoV-MA15-AE, or rSARS-CoV- 
MAIS. At 8 hpi, the culture medium was harvested and replaced with 
fresh medium supplemented with 2% fetal bovine serum (FBS). Culture 
supernatants were collected at 11 hpi (“nascent virus”), and virus titer was 
determined on Vero E6 cells. Viral RNA was isolated from 400 pl of 
supernatants using a MagMAX viral RNA isolation kit (Life Technolo- 
gies) according to the manufacturer’s instructions, and genomic RNA was 
quantified by reverse transcription-quantitative PCR (RT-qPCR) analy- 
sis. CDNAs were generated using a high-capacity cDNA transcription kit 
(Applied Biosystems) and the reverse primer Q-SARS-2015-RS (5’-ATG 
GCGTCGACAAGACGTAAT-3’). The cDNAs were amplified by qPCR 
using SYBR green PCR master mix (Applied Biosystems) with forward 
primer Q-SARS-1931-VS (5’-ACCACTCAATTCCTGATTTGCA-3’) 
and Q-SARS-2015-RS. The ratio of infectious virus titer to genomic RNA 
represents the specific infectivity. 

Cycloheximide chase to estimate the half-life of rSARS-CoV- 
MAI15-E* mutants. Subconfluent monolayers (90% confluence) of Vero 
E6 cells were infected with the wt or rsARS-CoV-MA15-E* virus at an 
MOI of 0.3. At 12 hpi, the cells were treated with 200 g/ml cyclohexi- 
mide (Sigma). Cells were harvested at 0, 2, 4, 6, and 8 h, and cell extracts 
were prepared as previously described (60). Samples were analyzed by 
sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) 
and Western blotting. In fact, E protein half-life measurements were ob- 
tained by calculating E protein relative mass normalized to actin mass for 
each time point. 

Western blot analysis. Cell lysates were resolved by SDS-PAGE, trans- 
ferred to a nitrocellulose membrane by wet immunotransfer, and pro- 
cessed for Western blotting. The blots were probed with a rabbit poly- 
clonal Ab specific for E protein (dilution of 1:6,000; kindly provided by 
Shen Shuo, Institute of Molecular and Cellular Biology, Singapore), a 
monoclonal Ab specific for hemagglutinin (HA) tag (dilution of 1:1,000; 
Sigma), or an antibody specific for B-actin (dilution of 1:10,000; Abcam). 
Bound antibodies were detected with horseradish peroxidase-conjugated 
goat anti-rabbit or anti-mouse antibodies (dilution of 1:30,000; Cappel) 
and the Immobilon Western chemiluminescent substrate (Millipore). 
The blots were quantitated by phosphorimaging. The half-life of each E 
protein mutant was calculated by plotting the signal intensity (ratio of 
band intensity). 

Plasmids. The pcDNA3-E plasmid encoding SARS-CoV E protein 
was used as previously described (25). To construct a pcDNA3 plasmid 
encoding the viral M protein fused to an HA tag (pcDNA3-HA-M), prim- 
ers M-EcoRI-HA-VS (5’-GCGCGCGCGAATTCGCCGCCATGTACCC 
ATACGATGTTCCAGATTACGCTGCAGACAACGGTACTATTACCG 
TTGAG-3’, encoding an HA sequence) and M-XholI-RS (5’-CGCGCTC 
GAGTTACTGTACTAGCAAAGCAATATTGTC-3’) were used. GAATTC 
and CTCGAG are the recognition sequences for EcoRI and Xhol, respec- 
tively. A DNA fragment encoding the SARS-CoV M protein was generated 
using as the template the plasmid pBAC-SARS-CoV-MAI15. To generate 
pcDNA3 plasmids encoding the E protein with deletions (pcDNA3-E- 
A3), primers A-EcoRI-VS (5'-GCGCGCGCGAATTCGCCGCCATGTA 
CTCATTCGTTTCGGAAGAAACAG-3’) and A-XhoI-RS (5’- 
CGCGCTCGAGTTAGACCAGAAGATCAGGAACTCCTTCAGAA 
GAGTT-3’) were used to amplify the DNA fragment encoding the SARS- 
CoV E deletion proteins using as the template plasmids pBAC-SARS- 
CoV-MA15-E* with the appropriate deletion. After PCR, these DNA 
fragments were digested by restriction enzymes EcoRI and Xhol and cloned 
into pcDNA3-E. All expression plasmids were verified by sequencing. 

Immunoprecipitation assay. Vero E6 cells were grown to 90% con- 
fluence and cotransfected with a plasmid encoding an N-terminal HA- 
tagged M protein (pcDNA3-HA-M), combined with plasmids expressing 
the full-length E protein or E protein with small deletions (pcDNA3-E or 
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FIG 5 rSARS-CoV-MA15-E* growth in the lungs of infected mice. BALB/c 
mice were intranasally inoculated with 1 X 10° PFU of the indicated viruses. At 
2 and 4 days p.i., lung tissue was harvested and viral titers were analyzed in 
Vero E6 cell monolayers. Means and standard deviations are shown (n = 3 
mice). 


pcDNA-E-A3). Cells cotransfected with the plasmid pcDNA3-HA-M and 
the empty plasmid pcDNA3, lacking the coding region of E protein, were 
used as controls. Twenty-four hours later, cell extracts were collected as 
previously described (60). For immunoprecipitation assays, monoclonal 
anti-HA agarose-conjugated clone HA-7 (Sigma) was used by following 
the manufacturer’s instructions. Briefly, 75 wl of the anti-HA agarose 
conjugate was incubated with the cell extracts overnight at 4°C. Immune 
complexes were eluted using 20 wl 2 SDS sample buffer and heating at 
95°C for 3 min. Analysis of precipitated complexes was carried out by 
SDS-PAGE and Western blotting. 

Virus infection and titration in mice. BALB/c mice were anesthetized 
with isoflurane and intranasally inoculated with 1 X 10° PFU of virus in 50 
wl of Dulbecco’s modified Eagle’s medium (DMEM). For protection ex- 
periments, mice were immunized intranasally with 6 < 10° PFU of the 
attenuated viruses and then challenged with an intranasal inoculation of 
1 X 10° PFU of rSARS-CoV-MAI5 at 21 days postimmunization. Mice 
were monitored daily for weight loss and mortality. SARS-CoV-MA15 
titers were measured in Vero E6 cells as described above. Viral titers were 
expressed as PFU/g tissue. 

Histopathologic examination of infected mouse lungs. Lungs were 
removed from infected mice, fixed in zinc formalin, and paraffin embed- 
ded. Sections were stained with hematoxylin and eosin. 

RNA analysis by RT-qPCR. BALB/c mice were intranasally inoculated 
with 1 X 10° PFU of the corresponding virus or DMEM, as a control. 
Lungs were removed and incubated in RNAlater (Ambion) at 4°C for 24h, 
prior to freezing at —80°C. Total RNA was extracted using an RNeasy 
minikit (Qiagen). Reverse transcription (RT) reactions were performed at 
37°C for 2 h using the high-capacity cDNA transcription kit (Applied 
Biosystems) to generate cDNAs. qPCRs were performed using TaqMan 
gene expression assays (Applied Biosystems) specific for murine genes 
(Table 2). Quantification was achieved using the 2 44CT method (61). The 
data represent the average of duplicate measurements from three inde- 
pendent mice. 

Microarray analysis. Lungs from infected mice were collected at 2 
days postinfection (p.i.), and total RNA was extracted using an RNeasy 
minikit (Qiagen). Total RNA from three different lungs was indepen- 
dently hybridized for each transcriptomic comparison. A total of 200 ng of 
RNA was amplified using a one-color low-input Quick Amp labeling kit 
(Agilent Technologies) and further purified with an RNeasy minikit 
(Qiagen). Preparation of probes and hybridization were performed as 
described in One-Color Microarray-Based Gene Expression Analysis, ver- 
sion 6.5 (62). Images were captured with an Agilent microarray scanner, 
and spots were quantified using feature extraction software (Agilent Tech- 
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FIG 6 Lung pathology caused by infection with rsARS-CoV-MA15-E* mutants. BALB/c mice were intranasally inoculated with 1 X 10° PFU of the indicated 
SARS-CoV deletion mutants and sacrificed at days 2 and 4 p.i. (A) Lungs were removed and sections were prepared and stained with hematoxylin and eosin. 
Asterisks indicate edema accumulation in both bronchiolar and alveolar airways. Original magnification was X20. Representative images are shown. (B) 
Macroscopic lung pathology in rsARS-CoV-MA15-E*-infected mice. Representative images of lung gross pathology. (C) Lung weight. Three lungs were 
evaluated in each case (n = 3 mice). Statistically significant data compared to mice infected with the wt virus are indicated (**, P < 0.01). 


2 DPI 


4 DPI 


4 
TIME POST-INFECTION, days 


3876 asm.org Journal of Virology April 2015 Volume 89 Number 7 


AYVHSEIT THH HSLSII avn Aq SLog ‘v2 Ainr uo /6so"wse'lAl//:djjy Wo. papeojuMOg 


nologies). Background correction and normalization of expression data 
were performed using linear models for microarray data (LIMMA) (63). 

Microarray data analysis. Linear model methods were used for deter- 
mining differentially expressed genes. Each probe was tested for changes 
in expression over replicates by using an empirical Bayes moderated t 
statistic (64). To control the false discovery rate (FDR), P values were 
corrected by using the method of Benjamini and Hochberg (64). The 
expected false discovery rate was controlled to be less than 5%. Genes were 
considered differentially expressed when the FDRs were <0.05 and fold 
changes were >2 or <-2. 

Measurement of leukocytes in the mouse lungs. Sixteen-week-old 
BALB/c mice were infected as indicated above with 1 X 10° PFU of each 
virus. Lungs were removed from mice at 4 days p.i., and cells were pre- 
pared as previously described (51). The numbers of macrophages, neu- 
trophils, and CD4 and CD8 T cells were determined as previously de- 
scribed (32, 58). 

Cytokine multiplex analysis. Lungs from infected mice were homog- 
enized and nuclear, cytoplasmic, and extracellular proteins were extracted 
using a nuclear extract kit (Active Motif, Carlsbad, CA). Supernatant pro- 
teins (cytoplasmic and extracellular proteins) were diluted 1:5 in assay 
diluent (Millipore) prior to analysis. The expression of tumor necrosis 
factor (TNF), MCP-1/CCL2, KC/CXCLI, IL-5, IL-6, IL-13, IFN-y, mac- 
rophage inflammatory protein la (MIP-1a)/CCL3, and MIP-1B/CCL4 
was determined as previously described (58). 

Statistical analysis. Student’s t test was used to analyze differences in 
mean values between groups. All results are expressed as means + stan- 
dard errors of the means (SEM). P values of <0.05 were considered sta- 
tistically significant. 

Microarray data accession number. MIAME-compliant results of the 
microarrays have been deposited in the Gene Expression Omnibus data- 
base (GEO [National Center for Biotechnology Information], accession 
code GSE59185). 


RESULTS 


Generation, growth, and characterization of rSARS-CoV- 
MA15-E* viruses. To identify the E protein regions responsible 
for rsARS-CoV-AE attenuation, a set of mutant viruses (rSARS- 
CoV-MA15-E*) were generated (Fig. 1A). In rSARS-CoV-MAI15- 
Mut 1, the E protein amino-terminal region was modified by in- 
troducing point mutations, instead of deleting amino acids, in 
order to maintain unaltered the primary sequence of the 3b pro- 
tein, as 3b and E genes partially overlap. Sequential or partially 
overlapping small deletions of 6 to 12 amino acids were intro- 
duced in the carboxy-terminal region of E protein to generate 
rSARS-CoV-MA15-A2, -A3, -A4, -A5, and -A6. The engineered 
viruses were rescued in Vero E6 cells, cloned by three rounds of 
plaque purification, and sequenced to confirm the presence of the 
desired mutations. Mutant E* and wt proteins were similarly lo- 
calized within infected cells as determined by immunofluores- 
cence by using antibodies specific for the E protein (Fig. 2 and data 
not shown). All the mutant E proteins partially colocalized with 
the ERGIC marker as previously observed for the full-length pro- 
tein (25). These data indicated that none of the regions deleted 
were essential for the subcellular localization of E protein and that 
E protein functions that were associated with its localization were 
not affected. SARS-CoV E protein promotes virus production 
since elimination of the E gene decreases virus growth from 10- to 
100-fold, depending on the cell type (28). To delimit E protein 
regions involved in SARS-CoV production, growth kinetics of 
rSARS-CoV-MA15-E* were analyzed in Vero E6 and Huh7.5.1- 
infected cells (Fig. 1B). Maximum viral titers were observed at 48 
to 72 hpiin Vero E6 cells and at 72 hpi in Huh7.5.1 cells, for all the 
viruses. Mut 1 and A6 viruses reached peak titers (>10’ PFU/ml) 
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FIG 7 Specific infectivity of rsSARS-CoV-MA15-E* mutants. Vero E6 cells 
were independently infected with each of the rsARS-CoV-MA15-AE, -A3, and 
-A5 and wt viruses, at an MOI of 0.3. At 8 hpi, the culture medium was 
harvested and replaced with fresh medium supplemented with 2% FBS. Cell 
supernatants were collected at 11 hpi (“nascent virus”), and the levels of 
genomic RNA and infectious virus titer were analyzed by RT-qPCR and plaque 
assay, respectively. The ratio of infectious virus titer (PFU/ml) to genomic 
RNA is depicted in the graph as percentage of specific infectivity. Means and 
standard deviations are shown (n = 3 mice). 


similar to those observed for rSARS-CoV-MA15 in both cell lines. 
In contrast, the A3 and A5 viruses grew to lower titers (around 10° 
and 10° PFU/ml in Vero E6 and Huh7.5.1 cells, respectively), sim- 
ilar to those observed for virus lacking full-length E protein (AE). 
The A4 virus grew to titers intermediate between titers of the wt 
and AE viruses, reaching titers close to 10’ and 10° PFU/ml in 
Vero E6 and Huh7.5.1 cells, respectively. The A2 virus showed the 
lowest titers (10° and 104 PFU/ml in Vero E6 and Huh7.5.1 cells, 
respectively). These data indicated that deletion of SARS-CoV E 
protein regions 2, 3, and 5 reduced virus growth to the same extent 
that the deletion of full-length E protein did, suggesting that these 
regions themselves were essential for E protein function. In a sub- 
sequent experiment, we studied the A3 and A5 viruses, because 
viruses grew to titers similar to those of rsARS-CoV-MA15-AE. 

Stability of rsARS-CoV-MA15-E* mutants in cell cultures. 
Prior to analyzing the virulence of rsARS-CoV-MA15-E*, their 
stability in cell culture was examined. The viruses were passaged 
daily 8 times in Vero E6 cells, followed by sequencing of the distal 
third of the genome, encompassing the S gene to the 3’ end. Only 
minor changes were detected in the viral sequences of one repre- 
sentative isolate, suggesting that these viruses were in general ge- 
netically stable, at least in tissue culture (Fig. 3). To analyze 
whether the mutations that the viruses incorporate after growing 
them in tissue culture cells are compensatory mutations, growth 
kinetics of the original and passaged viruses were compared by 
infecting Vero E6 and Huh7.5.1 cells and analyzing virus titers at 
different times postinfection. No significant differences in the 
production of infectious viruses were observed among the pas- 
saged viruses and the original viruses (data not shown), suggesting 
that the mutations that the viruses incorporate after growing them 
in cell culture are not compensatory mutations. 

Virulence of rsARS-CoV-MA15-E* mutants. It was previ- 
ously shown that a SARS-CoV lacking the full-length E protein 
was fully attenuated in vivo (28, 29, 32, 58). To evaluate the rele- 
vance of E protein regions in virus virulence, BALB/c mice were 
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FIG 8 E protein stability of rSARS-CoV-MA15* mutants. Stability of E pro- 
tein deletion mutants was evaluated in relation to that of the full-length E 
protein after infection of Vero E6 cells with rSARS-CoV-MA15-E* mutants. At 
12 hpi, the cells were treated with cycloheximide. Cells were harvested at the 
indicated times, and the amount of E protein mutants was determined by 
Western blotting. (A) Membranes were probed with anti-E protein and with 
R-actin-specific antibodies as a loading control. (B) The graph represents the 
values obtained after densitometry analysis. The percentage of protein remain- 
ing after cycloheximide addition is represented. Bars represent standard devi- 
ations of the mean (n = 3 mice). 


mock infected or infected with rsARS-CoV-MA15-E*, rSARS- 
CoV-MA15-AE, and rSARS-CoV-MA15. Clinical disease and sur- 
vival were evaluated for 10 days. Mice infected with the parental 
virus showed disease symptoms from 2 days p.i., reflected by leth- 
argy and ruffled fur. These mice rapidly lost weight, and 80% of 
them died by day 6 (Fig. 4). Animals infected with A6 virus, car- 
rying a deletion in the last segment of the E protein, lost weight 
and died in a fashion similar to those inoculated with the wild- 
type virus (Fig. 4). In contrast, mice infected with A2, A3, A5, and 
AE viruses showed no signs of disease (data not shown), did not 
significantly lose weight, and uniformly survived (Fig. 4). Mice 
infected with the Mut 1 and A4 viruses lost 15% and 10% of their 
initial weight, respectively, but recovered from 3 days p.i., and all 
of them survived (Fig. 4). These data showed that E protein re- 
gions 1 and 4 and especially 2, 3, and 5 were required to maintain 
the optimal replication or virulence of SARS-CoV. 
rSARS-CoV-MA15-E* growth in vivo. To evaluate whether 
any of the E protein regions were required for optimal virus 


3878 jvi.asm.org 


Journal of Virology 


INPUT IPaHA 

oO o 

4 7 
Ww Ww Ww Ww 
+ + + + + + 
= ¢ F = ¢ &¢€ 
<q << + <<< + 
x= x= = Ww = x= x= Ww 


FIG 9 Interaction of SARS-CoV E protein deletion mutants with M protein. 
Vero E6 cells were cotransfected with a plasmid pcDNA3 encoding the N-ter- 
minal HA-tagged M protein, combined with plasmids expressing the E-wt 
protein or the E protein with a small deletion, E-A3. Cotransfections of a 
plasmid encoding HA-M protein and ofa plasmid without E protein were used 
as controls. Cells were lysed and analyzed by Western blotting with specific 
antibodies for the E protein and HA (left) or subjected to immunoprecipita- 
tion with the monoclonal HA-specific antibody. The presence of E and M 
proteins was analyzed in the precipitated fractions using E- or HA-specific 
antibodies (right). 


growth in vivo, BALB/c mice were infected with rSARS-CoV- 
MA15-E*, rSsARS-CoV-MA15-AE, or rSARS-CoV-MA15, and vi- 
ral titers in lungs were determined at 2 and 4 days p.i. (Fig. 5). Peak 
titers were reached by all viruses at 2 days p.i. and decreased 
around 5- to 20-fold at 4 days p.i. in all cases. The parental virus 
showed the highest titers, 8 X 10” PFU/g of lung tissue at day 2 and 
5 X 10’ PFU/g at day 4. Titers were reduced when E protein was 
fully eliminated (10° PFU/g at day 2 and 5 X 10* PFU/g at day 4), 
confirming that rsARS-CoV-MA15-AE replicated less robustly in 
mice. rSsARS-CoV-MA15-A6 grew to 7 X 10” PFU/g at 2 days p.i. 
and 10’ PFU/g at 4 days p.i., similarly to the wt virus (Fig. 5), 
indicating that E protein carboxy terminus region 6 is the least 
relevant for optimal in vivo growth. The A2 and A5 viruses grew to 
titers similar to those observed for the AE virus, or even lower. In 
contrast, the Mut 1, A3, and A4 viruses grew to intermediate titers 
(10’ PFU/g and 10° PFU/g at days 2 and 4, respectively) in lungs. 
These results are similar to those observed in tissue culture cells 
and showed that especially regions 2 and 5 of SARS-CoV E protein 
are required for optimal virus growth. 

Lung pathology in rsARS-CoV-MA15-E*-infected mice. To 
analyze the mechanisms leading to rSsARS-CoV-MA15-E* atten- 
uation, macroscopic and microscopic changes in the lungs of in- 
fected mice were analyzed at 2 and 4 days p.i. (Fig. 6). The lungs 
from animals infected with the wt or A6 viruses but not AE, Mut 1, 
A2, A3, A4, and A5 viruses were increased in volume and were 
dusky red. Profuse areas of hemorrhage and weight increase, pos- 
sibly due to leukocyte infiltration and edema accumulation, were 
especially evident at 4 days p.i. (Fig. 6B and C). In contrast, histo- 
logical examination of lungs from AE, Mut 1, A2, A3, A4, and A5 
virus-infected mice revealed only a small amount of lung damage and 
cellular infiltration, whereas examination of lungs infected with wt or 
A6 virus showed substantial lung damage characterized by peribron- 
chial/perivascular inflammatory cell infiltration, edema, and thicken- 
ing of the alveolar walls at 2 and 4 days p.i. (Fig. 6A). These data 
demonstrated that E protein regions 1, 2, 3, 4, and 5 contribute to the 
lung damage induced after SARS-CoV infection. 

Specific infectivity and half-life of rSARS-CoV-MA15-E* 
mutants. In order to investigate whether the different small dele- 
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FIG 10 Leukocyte infiltrates present in the lungs of rsARS-CoV-MA15-E*-infected mice. BALB/c mice were intranasally infected with 1 x 10° PFU of the 
indicated SARS-CoV deletion mutants and sacrificed at 4 days p.i. The total number of leukocytes (A) and total numbers and percentages of macrophages (B), 
neutrophils (C), CD4* T cells (D), and CD8* T cells (E) were determined. Error bars represent the standard deviations of the means. Data are representative of 
three independent experiments (n = 3 mice). Statistically significant differences compared to infection with the wt virus are indicated (*, P < 0.05). 


tions within the E protein altered virus-specific infectivity, the 
titers of extracellular rsARS-CoV-MA15-wt, -AE, -A3, and -A5 
viruses were compared to the level of the corresponding extracel- 
lular genomic RNA. The specific infectivities of rSARS-CoV- 
MAI5-AE, -A3, and -A5 mutants were lower (21%, 11%, and 
17%, respectively) than that of the wt virus (Fig. 7). Thus, intro- 
ducing small deletions or modifications within E protein de- 
creases the infectivity of the particles, probably contributing to the 
lower titers observed for the rSARS-CoV-MA15-AE, -A3, and -A5 
mutants. 

To analyze whether the small deletions introduced in the E 
protein changed its stability, cellular translation was inhibited 
with cycloheximide, and the degradation of E protein was ana- 
lyzed by Western blotting. The half-lives of the E proteins encoded 
by the rSARS-CoV-MA15-A3 and -A5 virus mutants were longer 
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(3.5 h and 3.6 h, respectively) than that of the wt E protein, which 
had a half-life of 2.4 h, in agreement with previous results (65) 
(Fig. 8). The fact that wild-type E protein has a shorter half-life can 
be explained, because this protein is more abundantly expressed at 
early times compared to mutant proteins. These results suggest 
that the introduction of small deletions within the E protein mod- 
estly alter the stability of the protein. 

Binding between E protein and M protein. The interactions 
between SARS-CoV E and M proteins play an essential role in the 
viral particle assembly (66-68). Based on this observation, and on 
the result showing that the rSARS-CoV-MA15-A3 and -A5 viruses 
grew to lower titers than those for wt virus in cell cultures, the 
interaction between the E-A3 deletion protein and the M protein 
was analyzed. To this end, an HA-tagged, full-length M protein 
and E-A3 deletion protein were coexpressed in Vero E6 cells by 


jviasm.org 3879 


AYVHEIT THH HSLSII GVN Aq SLog ‘v2 Ainr uo /6so"wse'lAl//:dj}y Wos papeojuMOg 


AE vs mock . 


5 7911131517 5 7 911131517 5 7911131517 5 7 911131517 
AE vs wt A3 vs wt A5 vs wt 


17 dit 


LOG EXPERIMENT 


5 7911131517 5 7911131517 5 7 911131517 
LOG CONTROL 


B Cc 


Positive regulation of cell proliferation 
Positive regulation of cell death 
Regulation of cytokine production 
Translation 


FOLD CHANGE 
GENE AE vs WT A3 vs WT A5 vs WT 


INFLAMMATION 


CXCLS5 (Ena-78/Lix) -1.9 
In utero embryonic development CCL4 (Mip-1b)* -2.4 
po ee 
Blood | devel: it = 
cs cctsihietay 28 
‘ NFKB2 3.2 
CX3CL1 
0 104 10? 103 104 105 TNF* 
CCL2 (Mcp-1)* 
CXCL9 (Mig) 
CXCL13 
‘ 1L10 
Regulation of cell adhesion XCL1 (II8ra)* 
Di-, tri-valent inorganic cation homeostasis CSF3 (GCSF) 
Defense response CXCL10 (Inp10)* 
Positive regulation of cell adhesion feng 
Inflammatory response IL4* 
Response to wounding 1L13* 
40 40" 102 IMMUNE RESPONSE 
IL5* +4.4 +1.4 +6.4 
1L12B* 3.9 -2.3 3.6 


IFN RESPONSE 
Acute inflammatory response 


Extracelular structure organization 
Response to cytokine stimulus 
Endocytosis 

Inflammatory response 

Defense response 

Response to wounding 

Immune response 


Cis 
c1Qc 
SERPING1 


10 107 107 10% 10* 105 10% 
TRANSCRIPTION 
P-VALUE 


*Indicate the number genes by group STATS oe as a 


FIG 11 Differential gene expression in rSARS-CoV-MA15-E*-infected mice. Total lung RNA was extracted from mice infected with the indicated SARS-CoV 
deletion mutants at 2 days p.i. (n = 3 mice). (A) Differential gene expression was measured using microarrays. Only genes with an FDR less than 0.05 were 
considered candidate genes. Points with a differential expression higher or lower than 2-fold are represented as darker or lighter dots, respectively. (B) Candidate 
genes that were up- and downregulated in AE, A3, and A5 viruses were grouped according to Gene Ontology terms. Numbers on the x axis indicate DAVID FDR 
values. (C) Genes differentially expressed in AE-, A3-, and A5-infected mice compared to in wt-infected mice were classified according to their main biological 
functions. Black and gray lettering is used to indicate up- and downregulated genes, respectively. Asterisks indicates genes whose expression was confirmed by 
RT-qPCR. The numbers indicate the fold change for each gene in AE-, A3-, and A5-infected mice compared to in wt-infected mice. For those genes detected with 
more than one probe, the value corresponding to the highest upregulation or downregulation is represented. 
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TABLE 3 Analysis of host gene expression using mouse microarrays® 


Total no. of 

genes No. of No. of 

differentially upregulated downregulated 
Comparison expressed genes genes 
WT vs mock 2,918 1,473 1,445 
AE mutant vs mock 957 594 363 
A3 mutant vs mock 1,466 1,004 462 
A5 mutant vs mock 1,362 955 409 
AE mutant vs WT 1,382 752 630 
A3 mutant vs WT 83 39 44 
A5 mutant vs WT 286 44] 385 


* Upregulated, fold change of >2, FDR of <0.05; downregulated, fold change of <—2, 
FDR of <0.05. 


transient transfection, and coimmunoprecipitation experiments 
were performed. The E-wt and E-A3 proteins were specifically 
coimmunoprecipitated together with the M protein, using an an- 
tibody specific for the HA tag fused to the M protein (Fig. 9). These 
results confirm that E and M proteins interacted in infected cells 
and indicated that the A3 deletion introduced within the E protein 
did not significantly alter the binding of this protein to the M 
protein. 

Immune cell responses in rSARS-CoV-MA15-E*-infected 
mice. SARS progression and disease outcomes are associated with 
a massive influx of inflammatory cells to the lungs and a poor T 
cell response (50, 51, 54, 69). Leukocyte infiltrates in the lungs of 
mice infected with virulent rsARS-CoV-MA15-wt or representa- 
tive attenuated mutants (rSARS-CoV-MA15-AE, -A3, and -A5) 
were quantified and characterized by flow cytometry. Total leuko- 
cyte numbers were slightly increased in lung cell infiltrates of mice 
infected with the attenuated mutants (AE, A3, and A5 viruses), 
compared to those infected with the parental virulent virus (Fig. 
10A). Higher total and relative levels of macrophages were ob- 
served in lung infiltrates of the mice infected with the attenuated 
viruses than in those infected with rsARS-CoV-MA15 (Fig. 10B). 
Interestingly, a decrease in the percentage and total number of 
neutrophils was observed in the mice infected with the attenuated 
viruses compared to those in rsARS-CoV-MA15-infected mice 
(Fig. 10C). These data showed that the attenuation observed for 
mutant viruses (AE, A3, and A5 viruses) correlated with a decrease 
in neutrophil influx and an increase in macrophage influx into the 
lungs. Higher total and relative numbers of CD4* and CD8* T 
cells were present in the lungs of AE, A3, and A5 virus-infected 
mice than in virulent rSARS-CoV-MA15-infected lungs (Fig. 10D 
and E). All data analyzed were statistically significant, with the 
exception of total leukocyte, neutrophil, and CD4* numbers. 
These data suggested that mice infected with the attenuated vi- 
ruses develop higher T cell responses than those infected with 
virulent SARS-CoV, most likely contributing to virus clearance 
and disease attenuation. 

Effect of SARS-CoV-MA15-E* virus infection on host gene 
expression. To analyze host responses responsible for the at- 
tenuation of the rsARS-CoV-MA15-E* mutants, the transcrip- 
tomes of rsARS-CoV-MA15-, rSsARS-CoV-MA15-AE-, rSARS- 
CoV-MA15-A3-, rSARS-CoV-MA15-A5-, and mock-infected 
lungs were compared, using mouse microarrays (Fig. 11A and 
Table 3). The genes differentially expressed during SARS-CoV 
infection (both up- and downregulated) in rsARS-CoV-MAI15- 
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AE, -A3, and -A5 virus-infected mice, compared to wt-infected 
mice, were clustered using DAVID software (70). The most statis- 
tically significant clusters were associated with the immune re- 
sponse, the response to wounding, the defense response, and the 
inflammatory response (Fig. 11B). Based on their principal bio- 
logical function, the genes downregulated in mice infected with 
the attenuated rsARS-CoV-MA15 mutants were involved mainly 
in inflammation (Fig. 11C). Interestingly, this reduction corre- 
lated with relative increase in the expression of the Th2 genes, such 
as those encoding IL-4 and IL-13, in comparison to that in mice 
infected with the parental virus. 

To validate and confirm the results obtained using microar- 
rays, several genes involved in inflammatory responses were se- 
lected for RT-qPCR analysis, and their expression was tested at 2 
and 4 days p.i. A clear reduction of all the proinflammatory cyto- 
kines analyzed (TNF, CCL2, CCL3, CCL4, CXCL1, CXCL2, 
CXCLI1O, IL-6, IL-12B, and IFN-y) was observed in the lungs of 
mice infected with the AE, A3, and A5 attenuated viruses, com- 
pared to that in wt-infected mouse lungs (Fig. 11C and 12A and 
B). In contrast, the expression of the Th2 cytokines IL-4, IL-5, and 
IL-13 was upregulated in the lungs of mice infected with the AE, 
A3, and A5 attenuated viruses, compared to those from rSARS- 
CoV-MA15-infected mice (Fig. 11C and 12A and B). All data an- 
alyzed were statistically significant with minor exceptions. These 
changes were specific, as the mRNA levels of rRNA 18S and trans- 
forming growth factor B (TGF-B), encoded by another gene in- 
volved in the anti-inflammatory response but regulated posttran- 
scriptionally (71, 72), were the same in all mice (Fig. 12B). 

To confirm that the mRNA levels correlated with protein ac- 
cumulation, the levels of TNF, CCL2, CXCLI1, IL-5, IL-6, IL-13, 
IFN-y, CCL3, and CCL4 were evaluated in lung extracts from 
mice infected with the rSsARS-CoV-MA15, A3, A5, and AE viruses 
or mock infected, as a control. mRNA levels of these cytokines 
matched with their protein levels in the mouse lungs (Fig. 13). All 
cytokine differences at the protein level were statistically signifi- 
cant, with minor exceptions in CCL3, CXCL1, and CCL4 cyto- 
kines. All together, these results indicated that the macroscopic 
and microscopic inflammation present in the lungs of infected 
mice correlated with overexpression of proinflammatory cyto- 
kines and downregulation of anti-inflammatory cytokines. 

Immunization with the attenuated mutant viruses protects 
mice against challenge with SARS-CoV-MAI5 virus. Live atten- 
uated vaccines are considered most effective in their ability to 
induce a long-lived, balanced immune response. In order to de- 
termine whether the attenuated mutant viruses were good vaccine 
candidates, the induction of protection conferred by the attenu- 
ated viruses with point mutations or small deletions, against the 
challenge with virulent SARS-CoV-MA15, was studied. Groups of 
BALB/c mice were either mock immunized with PBS or immu- 
nized with the Mut 1, A2, A3, A4, and A5 attenuated viruses. At 21 
days postimmunization, mice were challenged with rsARS-CoV- 
MAI15, and weight loss and mortality were monitored for 14 
days postchallenge (Fig. 14). All the nonimmunized mice lost 
weight and died by day 7 after rSARS-CoV-MA15 challenge. In 
contrast, vaccination with the attenuated mutant viruses com- 
pletely protected mice from the lethal dose of SARS-CoV- 
MAI1S5, as they showed no weight loss and all survived. These 
data indicated that all attenuated mutant viruses effectively 
induced a high level of protection against lethal rSARS-CoV- 
MA15 challenge in mice. 
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compared to mice infected with wt virus are indicated (*, P < 0.05). 


DISCUSSION 


An extensive study using a collection of recombinant viruses was 
performed to analyze E protein regions and host responses that 
contributed to the virulence of SARS-CoV. Several mutants en- 
coding small deletions or point mutations within E protein were 
engineered, in the context of mouse-adapted SARS-CoV-MA15 
(57). Point mutations in the amino-terminal region and small 
deletions in the internal carboxy-terminal region of E protein led 
to viruses with an attenuated phenotype, similarly to that obtained 
by deleting the whole E protein of SARS-CoV. Infection of mice 
with the rSsARS-CoV-MA15-Mut 1, -A2, -A3, -A4, -A5, and -AE 
mutants led to decreased lung pathology compared to that of wt- 
infected mice and no mortality. Therefore, these regions are re- 
quired for the virus to develop its full pathogenic potential. In 
contrast, virus lacking the 12 carboxy-terminal amino acids of E 
protein (A6 mutant) was fully virulent, reinforcing the idea that 
not all small deletions within the E protein were attenuating. The 
deletion within rsARS-CoV-MA15-A6 includes the PDZ-binding 
motif (PBM) (73, 74), and we previously showed that a mutant 
lacking this region is attenuated (73). One possible explanation for 
this apparent discrepancy is that the last 6 amino acids within the 
A6 virus (YSRVKN; Fig. 1) form an alternative PBM that could 
restore the virulent phenotype of the virus, an aspect that is cur- 
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rently being explored in our laboratory. Deletion or mutation of E 
protein regions could in principle affect basic properties of the 
protein, such as its subcellular localization and half-life, which 
could broadly disturb its major functions within the cell. The cy- 
toplasmic tail of SARS-CoV E protein was shown to contain a 
signal providing Golgi localization (75). In addition, CoV E pro- 
tein may impair protein trafficking depending of its half-life (65). 
None of these aspects was significantly altered, as mutant E pro- 
teins colocalized mainly with ERGIC markers, as described for the 
wt protein (25), and the half-life of E* protein mutants was slightly 
increased (1.5-fold) compared with that of the wt protein (65). 
Furthermore, the interaction between E and M proteins necessary 
for virus morphogenesis (66-68) was preserved despite deleting 
the A3 region. Thus, the E-A3 region of SARS-CoV E protein is 
not required to interact with M protein. In contrast, the specific 
infectivity of the different small deletions within the E protein was 
lower (6-fold) than that of the wt virus. These modifications in E 
protein possibly contributed to the reduced infectivity observed 
for E protein deletion mutants that exhibited an attenuated phe- 
notype. However, these phenomena are not necessarily linked, 
because we previously demonstrated that recombinant viruses 
containing point mutations in the E protein that inhibited its ion 
channel activity or eliminated its PBM grew to titers similar to 
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Weight loss (A) and survival (B) were recorded daily. 


those of the parental virus in the lungs of the infected mice (37, 
73). Nevertheless, these viruses caused a reduction in lung inflam- 
mation and pathology and were attenuated (37, 73). Therefore, E* 
protein mutants may decrease virus virulence by affecting specific 
host responses, independently of effects on virus replication. In- 
terestingly, the patterns of gene expression in infected lungs re- 
vealed that proinflammatory genes, such as those encoding TNF, 
CCL2, CCL3, CCL4, CXCL1, CXCL2, CXCL10, IL-6, IL-12B, and 
IFN-y, were downregulated in rsARS-CoV-MA15-AE, -A3, and 
-A5-infected mice compared to in wt-infected mice, most proba- 
bly contributing to viral attenuation. Notably, increased levels of 
proinflammatory cytokines such as TNF, IL-6, CXCL10, CCL2, 
IL-8, and IL-12, which play a major role in mediating and ampli- 
fying ALI/ARDS, were found in the blood and the lungs of SARS 
patients (40-45, 48, 76). Furthermore, the expression of the Th2 
cytokines IL-4, IL-5, and IL-13 was increased in the lungs of mice 
infected with the attenuated viruses compared to in those infected 
with the virulent wt virus, possibly contributing to a better out- 
come. In agreement with this finding, the expression of IL-13 was 
increased in asymptomatic young BALB/c mice infected with a 
mouse-adapted SARS-CoV, compared to aged infected mice with 
severe respiratory disease (77). All these data strongly suggest that 
the lower proinflammatory cytokine expression and higher anti- 
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inflammatory cytokine levels observed after infection with the at- 
tenuated mutants could account for the diminished lung pathol- 
ogy and virus attenuation. 

Neutrophils play a key role in the progression of ALI/ARDS, 
and SARS patients with adverse outcomes presented with neutro- 
philia and high neutrophil counts in the lungs (38, 54, 78, 79). 
Similarly, massive neutrophil influx into the lungs contributed to 
disease after influenza A virus and respiratory syncytial virus 
(RSV) infections (80-83). This phenomenon was reproduced in 
mice infected with virulent SARS-CoV-MA15, whereas a signifi- 
cant reduction of neutrophils in the lungs of the mice infected 
with the attenuated mutants was observed. Accordingly, the ex- 
pression of CXCL1 and CXCL2, key factors that induce the che- 
motaxis of neutrophils (84, 85), was increased when the full- 
length E protein was present during SARS-CoV infection. 

An increase of macrophage numbers was observed in mice 
infected with the attenuated mutants, compared to those infected 
with the virulent wt virus. In addition to their role in immunity, 
macrophages maintain physiologic homeostasis and contribute to 
tissue repair (86, 87). Therefore, an increase of these cells during 
infection could favor disease resolution. Together, these data in- 
dicate that an increase in neutrophil numbers in conjunction with 
a decrease in macrophage numbers in the lungs is associated with 
more severe SARS-CoV pathogenesis and that the deletion of E 
protein carboxy-terminal regions was responsible for these 
changes in the phenotype of infiltrating immune cells and virus 
attenuation during SARS-CoV infection. 

T cells are important for viral clearance and for protection 
from disease in primary SARS-CoV infection (50-52), and low T 
cell counts correlate with adverse outcomes in SARS patients (38, 
53, 54, 56). Interestingly, we have shown that in rsARS-CoV- 
MA15-AE, -A3, and -A5-infected mice, which do not show a sig- 
nificant disease, CD4* and CD8* T cell counts are higher than 
those found in mice infected with the parental virus, which show 
adverse outcomes after infection. We were unable to measure vi- 
rus-specific T cell responses, because rsARS-CoV-infected mice 
die before a significant T cell response develops (Fig. 4). Several 
anti-SARS-CoV vaccine candidates have been developed since 
2003, but none has been approved for use in humans (88, 89). Live 
attenuated vaccines are considered most effective compared to 
other types of vaccines, as they induce a long-lived, balanced im- 
mune response. Remarkably, the attenuated mutants generated in 
this work provided full protection after challenge with the virulent 
wt virus. These viruses express a complete repertoire of viral pro- 
teins, except for the small regions deleted or mutated in the E 
protein. Therefore, it is predicted that they will elicit both anti- 
body and T cell responses directed against the virus. The major 
problems of using live attenuated vaccines are the possibility that 
the viruses may revert to virulence. To overcome this limitation, 
we are introducing additional attenuating mutations into the 
Nsp1 protein to generate a safer vaccine candidate. The Nsp1 gene 
was selected as a target because it is located at a very distant site 
from that of the E gene in the viral genome, making the generation 
of a virulent virus through single recombination events with cir- 
culating coronaviruses unlikely. 

In summary, we have shown that changes in the amino termi- 
nal and internal carboxy-terminal regions of E protein led to sta- 
ble attenuated SARS-CoV mutants. These mutants showed a 
slight reduction in their specific infectivity, maintained the same E 
protein subcellular localization, and slightly increased their half- 
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life compared to the wild-type virus. Moreover, the A3 deletion 
introduced within the E protein maintained the interaction be- 
tween E and M proteins. SARS-CoV attenuation correlated with 
limited lung pathology, mediated by decreased proinflammatory 
cytokine expression, increased anti-inflammatory cytokine levels, 
and reduced neutrophil accumulation in the lungs. In addition, 
higher CD4* and CD8”* T cell counts in mice infected with the 
attenuated mutants most likely contributed to better outcomes. 
These attenuated mutants fully protected mice from challenge 
with virulent virus. These studies are the basis for the rational 
design of SARS-CoV vaccine candidates. 
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